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Abstract 

This work provides a quaternionic representation for real symplectic matrices in dimension 
four, analogous to that for the orthogonal group. It also provides a technique to compute 
this representation from the entries of the matrix being thus represented. As a byproduct it 
shows how one can compute the polar decomposition of a 4 x 4 symplectic matrix without 
any recourse to spectral calculations. 

PACS Numbers: 03.65.Fd, 02.10.Yn, 02.10.Hh 



1 Introduction 

The most important groups in physics are, arguably, the rotation, and the symplectic groups, 
[1,6]. In particular, the symplectic group is central to classical mechanics, [14], to classical and 
quantum optics, [1], quantum mechanics and quantum information processing, [3]. Therefore, 
having as many parametrizations of these groups as possible is desirable. 

In dimension four, there is a well known parametrization of the orthogonal group via 
a pair of unit quaternions. This has innumerable applications in physics and engineering, 
[13, 15]. For the symplectic group, S'p(4, R), there seems to be none. Note, by the "symplectic 
group" , we refer to the real symplectic group, and not the similarly named group Sp{n), which 
preserves the standard inner product on i/" (here H stands for the quaternions). The latter 
is, of course, defined via quaternions, [15]. It is the intention of this note to provide such a 
parametrization for 5p(4, R). Since the symplectic group is 10 dimensional (in comparison to 
the six dimensions of the 30(4, R)), the parametrization to be provided is accordingly more 
intricate. 

For motivating the balance of the paper (in particular, the contents of the introductory 
section) , we assume temporarily some familiarity with the algebra ismorphism between H(i)H 
(the tensor product of the quaternions, H, with itself) and the algebra of 4 x 4 real matrices, 
denoted Mi{R), [15, 12, 5, 9, 16, 17, 19, 20]. See section 2 for more details. Using this 
algebra isomorphism, it would seem that to obtain a quaternion representation of an X G 
Sp{A, R), a logical approach would be to write out the conditions imposed on the quaternion 
representation of an X G Mi{R) by the relation 

where J4 is the defining matrix of the symplectic group (its quaternion representation is 10 j). 
However, this produces an immensely complicated system of equations. 

Indeed, if one were to pursue this approach for the orthogonal group, i.e., write out the 
condition X^X = 7 in quaternion form, one does not recover the pair of unit quaternions 
representation. This is not surpirising since those orthogonal matrices with determinant equal 
to —1 are not represented by decomposable tensors, i.e., by an element of H®H, of the special 
form u^v, with u,v G H. Only if one were to insist that the quaternion representation of the 



X one began with was to be already given by a decomposable tensor, would one recover the 
fact that an X e <S'0(4, R) is represented by a tensor u®v, with u and v unit quaternions. Note 
that the condition, Aet{X) = 1 for the group SO{A,R) is subsumed by this representation, 
since the determinant of a matrix represented by p (8> ?, with p,q G H, is 

ibinuir 

as a calculation shows. 

However, it is not immediately evident that a special orthogonal matrix must be repre- 
sented by a decomposable tensor. In our opinion, the simplest way to see this is to use the fact 
that the exponential map from so(4, R) (the Lie algebra of 4 x 4 anti-symmetric matrices) is 
onto 50(4, R). Since the quaternion reprsentation of anti-symmetric matrices is easy to find 
(see section 2), a calculation, [19], shows that indeed matrices in 50(4, R) must be reprsented 
by decomposable matrices. Unfortunately, the analogous statement, viz., that the exponen- 
tial map from sp(4, i?) (the Lie algebra of Hamiltonian matrices) is onto Sp{4:,R), is false. 
Even if it were to be true, the exponential of Hamiltonian matrices do not typically admit a 
decomposable quaternion representation. Thus one cannot make the analogous simplification 
for the symplectic group. 

To obtain a different perspective on the complexity of the question at hand, we note 
that it is straightforward to get quaternionic representations of Hamiltonian matrices. This 
is because: i) Hamiltonian matrices are expressible as J2nS, with S real symmetric, and ii) 
quaternionic representations of symmetric matrices, and J4 are easy to find. Now Hamiltonian 
matrices form the Lie algebra of symplectic matrices. So it is reasonable to inquire if a 
similar translation from orthogonal to symplectic matrices is valid. However, this is patently 
incorrect, i.e., it is not true that a symplectic matrix can be represented as J4G, with G 
orthogonal. There also seems to be no easy modification of this approach which will do the 
job. 

To understand, the approach taken in this work, we note that the initial suggestion that 
one write out the quaternionic version of the condition X'^J4X = J4 leads to something 
tractable, provided one imposes the additional restriction that the original X be also real 
symmetric. This is very useful since it is known that in the polar decomposition X = UP, 



of a symplectic X, both the orthogonal factor (f/) and the positive definite factor (P) are 
symplectic as well (see [6, 14, 18], for instance). Since P is positive definite it is symmetric, and 
this explains the utility of starting with the assumption that X is symmetric and symplectic. 

For a variety of reasons, it is more convenient to obtain quaternionic representations of 
symmetric and symplectic X and then further refine these to obtain quaternionic character- 
izations of positive definite and symplectic X. Therefore, we eschew one route to the latter, 
viz., using the fact that a positive definite symplectic X must be the exponential of a sym- 
metric, Hamiltonian matrix (this exponentiation can be performed in closed form, [19] - but 
it would not yield the maximal benefit). Characterizing symmplectic, orthogonal matrices 
is relatively easy. Combining the two one gets a quaternionic representation of symplectic 
matrices. These are stated in Theorem (3.1) and Theorem (3.4), with the latter being an 
improvement of the former which takes into account the characeterization of positive definite, 
symplectic matrices obtained in Theorem (3.3) . 

Let us now consider the question of finding formulae for these quaternions in terms of 
the entries of the given symplectic matrix X. The first step is to find P. Now P is not any 
arbitrary symplectic, positive definite matrix. It is, as in the case of all polar decompositions, 
the unique positive definite square root of X^X. The traditional route to finding P involves 
diagonalizing X'^ X [see Remark (3.2)]. In this work we avoid this. Rather, we provide an 
explicit algorithm for finding P which requires only the solution of a very simple 2x2 linear 
system of equations. This is achieved by explicitly showing that there are only two real sym- 
metric, symplectic matrices H, with positive trace, satisfying H"^ = X^X. Since one of these 
must be P, we can easily find P. For all this it is crucial to characterize symmetric (not neces- 
sarily positive definite) symplectic matrices via quaternions, which explains the route chosen 
in this work. Of course, these results on symmetric, symplectic matrices are of independent 
interest. We note, however, that this work does indeed obtain a quaternionic characterization 
of positive-definite, symplectic matrices [see Theorem (3.4)]. In particular, a quaternionic 
parametrization of two-mode squeezing operators is thereby obtained, since positive-definite 
symplectic matrices are precisely those matrices which represent such operations, [1]. 

There are, of course, other global factorizations of the symplectic group, such as the Euler 
(Cartan) and Iwasawa decompositions, [1, 6, 7, 8, 21], which could have been used as starting 



points for finding a quaternionic representation. But we found the polar decomposition as 
the most useful, since the polar decomposition of a matrix has innumerable applications, [10]. 
Thus, the fact that the polar decomposition based representation in this work, also yields 
a simple, constructive procedure for finding the polar decomposition itself as a byproduct 
played an important motivating role. 

It is appropriate at this point to record some history of the linear algebraic applications 
of the isomorphism between H ® H and M4,{R). This isomorphism is central to the theory 
of Clifford algebras, [15]. However, it is only relatively recently been put into use for linear 
algebraic (especially numerical linear algebraic) purposes. To the best of our knowledge 
the first instance seems to be the work of [12], where it was used in the study of linear 
maps preserving the Ky-Fan norm. Then in [9], this connection was used to obtain the 
Schur canonical form explicitly for real 4x4 skew-symmetric matrices. Next, is the work of 
[5, 16, 17], wherein this connection was put to innovative use for solving eigenproblems of 
a variety of structured 4x4 matrices (including symmetric matrices). Finally, in [19, 20], 
this isomorphism was used to explicitly calculate the exponentials of a wide variety of 4 x 4 
matrices. 

The balance of this manuscript is organized as follows. In the next section some notation 
and preliminary results on symplectic matrices, positive definite matrices and the algebra iso- 
morphism between H^H and Mi{R) are collected. The next section contains all the results 
in this work. Proposition (3.1) provides the quaternion characterization of symmetric and 
symplectic matrices, central to this work. This is used in Theorem (3.1) to provide a quater- 
nion representation of Sp{A,R), which is further refined in Theorem (3.4). This refinement 
is based on a quaternionic characterization of positive definite matrices in Theorem (3.3). 
The principal technical tool in the proof of Theorem (3.3) is that of squaring a symmetric, 
symplectic matrix. This is also a key ingredient in the proof of Theorem (3.2). This latter 
theorem provides an explicit technique to calculate the polar decomposition of matrices in 
5^(4, R). While the proof of this theorem is somewhat lengthy, the key point is that it yields 
an algorithm for finding the positive definite factor in the polar decomposition, which requires 
only the calculation of a few inner products in R^ and the solution of a 2 x 2 linear system. 
This is summarized in an algorithm. Next a different perspective is provided on the paucity 



of symplectic, symmetric square roots of X^X, for an X G Sp{4:,R). As a byproduct an 
explicit formula for the characteristic polynomial of an X G 5*^(4, R) is obtained. The final 
section offers some conclusions. 



2 Notation and Preliminary Observations 

The following definitions, notations and results will be frequently met in this work: 
• M4{R) (also denoted gl{4:, R)) is the algebra of real 4x4 matrices. 



• J2n is the 2n x 2n matrix which, in block form, is given by J2n 



I , \ 



Sp{2n,R) denotes the Lie group of symplectic matrices, i.e., those 2n x 2n matrices, 
satsifying X^J2„X = J2n- sp{2n,R) is its Lie algebra, consisting of those real (resp. 
complex) 2n x 2n matrices which satisfy X'^J2n + J2n^ = 0. Such matrices are also 
called Hamiltonian. 

• If X e Sp{2n,R), then X'^ = -J2nX^J2n- Furthermore, if X G Sp{2n,R) then X^ is 
also in Sp{2n, R). 

• Essential use of the following theorem will be made in this work (see [6, 14, 18]): 

Proposition 2.1 Let X be a real symplectic matrix, and let X = PQ be its polar 
decomposition, with P positive definite and Q real orthogonal. Then P and Q are also 

real symplectic. 

See [18] for examples of other matrix groups for which an analogous statement holds. It 
is worth recalling here that P is the unique positive definite square root of the positive 
definite matrix X^X. 

• For a polynomial P{x) = X^"=o degree at most n, its reverse is the poly- 
nomial Pj:ev{x) = Yl^=o ^n-iX'^ ■ In this work we will use the fact that the charac- 
teristic polynomial of a symplectic matrix equals its reverse [see [18], for instance]. 
This follows from the well-known fact that for any invertible matrix m x m A, with 
characteristic polynomial p{x), the characteristic polynomial of is the polynomial 



(— l)™(det(A))'"a;'"p(i). Now for any polynomial, x'^p{^) is precisely its reverse. Fur- 
ther for symplectic A, m is even and det{A) is 1. Finally, a symplectic A is similar to its 
transpose. Hence, the characteristic polynomial of A is the same as that of its inverse, 
whence it equals its reverse. 

We next collect some definitions and results on real positive definite matrices. All details, 
together with extensions to the complex positive semidefinite case, may be found in [10]. 

• Definition 2.1 Let y be a real positive definite matrix. A real square matrix Z satis- 
fying Y = Z is said to said to be a square root of Y . 

See [10] for details (extending to the case of complex positive semidefinite matrices). 
Square roots of positive definite matrices are not unique. However, if Z\ is a square 
root of Y then Zi is also a square root of Y iff there exists a real orthogonal matrix U 
such that Zi = Z\U . 

• Let y be a real positive definite matrix. Then there exist real symmetric matrices H 
such that y = B."^ . Clearly any such if is a square root in the sense of Definition (2.1). 

• Let y be a real positive definite matrix. Then there exists a unique real positive definite 
matrix P with Y = . Clearly this P is an example of a real symmetric matrix whose 
square equals Y . 

Next relevant definitions and results regarding quaternions and their connection to real ma- 
trices will be presented. Throughout H will be denote the skew-field ( the division algebra) of 
the quaternions, while P stands for the purely imaginary quaternions, tacitly identified with 

H H and M4{R}: The algebra isomorphism between H ® H and 5^(4, R), which is central 
to this work is the following: 

• Associate to each product tensor q ^ H ® H, the matrix, Mp^q, of the map which 
sends x G H to pxq, identifying with H via the basis {l,i,j,k}. Thus, if p = 



Po + Pli + P2j + Psk; q = qo + qii + 92i + 93^, then 



Uo Vo Wo Zo 



Ui Vi Wi Zi 



U2 V2 W2 Z2 



\ 



Us V3 W3 Z3 



I 



with 



pq 



I + Uyi + U23 + Uzk 



piq = 



Vo + vii + V23 + v^k 



pjq = 



Wo + wii + W2j + wsk 



pkq 



- + zii + Z2j + zsk 



Here, q = qo - qii - q^j - qsk 

Extend this to the full tensor product by linearity, e.g., the matrix associated to 2(pi (g) 
qi) — 9{p2 <8) Q'2) is the matrix 2Mpj0gj — dMp^^q^- This yields an algebra isomorphism 
between H ®H and Mi{R). 

In particular, a basis for gl{A,R) is provided by the sixteen matrices Me^^ey as ex,ey 
run through l,i,j,k. Of these only M101 is not traceless. In particular, J4 = Mi^j 
belongs to this basis. 

Define conjugation in H ^ H hy first defining the conjugate of a decomposable tensor 

a eg) 6 as a (g) 6, and then extending this to all of ff (g) ff by linearity. Furthermore 



Thus, the most general element of M4[R) admits the quaternion representation al (g 
l+p(8>i + g(8>i + r(8>fc + s(8)l + l(8>t, with a £ R and p, q, r,s,t £ P. The summand 
al^l+p^i + q<S>j + r<Sik is the symmetric part of the matrix, while the summand 
s(g)l + l(8)iis the skew-symmetric part of the matrix. Expressions for a,p,q,r,s,t 
(which are linear in the entries of the matrix being represented) are easy to find, [16]. 
Finally 4a is the trace of the matrix. 



3 Quaternion Representations of 5^(4, R) 



In this section the main results of this work are presented. 

To develop a quaternionic representation of an X e 5^(4, R), let us invoke Theorem (2.1). 
Per that result it sufSces to obtain quaternionic representations of positive definite, symplectic 
matrices and symplectic, orthogonal matrices. 

Let X = UP he the polar decomposition of X. Since U is symplectic and orthogonal 
it must, in fact, be special orthogonal. Obtaining the quaternionic representation of such a 
matrix is easy. It is given hy q = u ^ v, with u, v unit quaternions (corresponding to special 
orthogonality) with the further restriction that vj = jv (corresponding to symplecticity) . 
This imposes no further restriction on u, but forces v to be of the form v = vq + V2j, with 
vl + v^ = l. 

To obtain a quaternionic representation of P, it is convenient (for the reasons mentioned 
in Section 1) to obtain a quaternionic characterization of symmetric, symplectic X's. This is 
achieved in: 

Proposition 3.1 Let X be a 4 x 4 symplectic matrix which is also symmetric. Then it 
admits the quaternion representation X = al<Sil+p<Sii + q®j + r<Sik, with aq = r x p, 
p.q = = r.q, and a satisfying the constraint — p.p + q.q — r.r = 1. If a, which is j'th the 
trace of X, is not zero then X is symplectic iS aq = r x p and — p.p + q.q — r.r = 1. 

Proof: The proof proceeds by equating the quaternion expansion of J4X = XJ^X, 
viz., (a l®l+p®i + q^j + r^k){l®j){a l®l+p®i + q®j + r®k) with 1 (S) j. 

This calculation is facilitated by the observation that for any matrix X (not necessarily 
symplectic or symmetric), the matrix X^J2„X is always skew-symmetric. Therefore, in the 
expansion (a 1 ® 1 + p ® i + q® j + r ® k){l ® j){a l®l-\-p®i + q®j-\-r®k) one needs to 
inspect only those terms of the form s 1 and 1 with s,t G P (i.e., those corresponding 
to the skew-symmetric terms) and equate all but one of them to zero. The non-zero term is, 
of course, the 1 <8) j term and this is equated to 1. Specifically, expansion of X^JiX is 

(a^ - p.p + q.q - r.r)(l ® j) + {2rxp- 2aq) ®1 + {2p.q)l ® i + (2r.g)l 18) k 

Note that the vanishing of the second term yields the condition aq = r x p. Thus, if a 9^ 0, 



then either q = or q = and this, of course, ensures the conditions p.q = = r.q. This 
finishes the proof. 

Now note that for the situation at hand, where X is not merely symmetric but, in fact, positive 
definite we necessarily have a ^ 0. This yields a preliminary result about the quaternion 
representation of symplectic matrices, which will be refined later in this section: 

Theorem 3.1 Let X be an element of Sp{4:,R). Then there exist a,vo,V2 £ R, p,q,r G P, 
and a unit quaternion u satisfying the constraints — p.p + q.q — r.r = 1, q = and 
' o + I'i = 1 such that X admits the following quaternion representation X = [u (Si {vq + 

V2j)] [a lig)l+p®i + gig)j+rig)fc]. 

Note that this result gives the correct dimension count for S'p(4, R). u is determined by 3 real 
quantities and Vo,V2 together yield one additional parameter. The non-compact portion of 
Sp{4:, R) is determined by the two vector quantities p, r, and the scalar a, with one restriction, 
thereby yielding 6 real parameters. Put together we get the dimension count to be 10. 

We now turn to the issue of computing these quaternions from the entries of X. Assuming the 
symmetric portion of the representation is available, the calculation of the factor [u'E){vo+V2j)] 
is amenable to the any technique which will yield the quaternion representation of a matrix 
in S0{4:,R). This ought to be folklore, but surprisingly the only explicit record of this that 
we were able to find is in [2]. 

Next, consider computing the symmetric, [a l(g)l+p0i-|- j + r A;] factor. This 
is, of course, not just any symmetric matrix. It is the unique positive definite square root 
of X'^X. To that end, one could, of course, resort to diagonalizing the symmetric matrix 
X^X. Instead, as we will see below that there is an even more explicit method for finding 
the symmetric factor, which reveals some facts which are of interest in their own right. 

Specifically, recall from Section 2, that any positive definite matrix is the square of a 
real symmetric matrix. The real symmetric matrix is not unique, but one of these is the 
unique positive definite square root of the positive definite matrix in question. Now, when 
the positive matrix is X^X^X £ S'p(4, i?), we know that the unique positive definite square 



root is also symplectic. Furthermore, being positive definite its trace is positive. So, instead 
of looking at all possible real symmetric matrices whose square is X^X, one needs to inspect 
only those which are symplectic and have positive trace, in addition. One then finds the 
pleasant conclusion that there are only few such candidates. 

Theorem 3.2 Let X e Sp{A,R). Then there are atmost two (and, at least one) matrices H 
which satisfy i) H'^ = X'^X; ii) H is real symmetric and symplectic; and iii) Trace {H) > 0. 
One of these is precisely the unique positive definite square root of X'^X, and thus the positive 
definite factor in the polar decomposition of X. Furthermore, H can be found explicitly via 
the solution of a simple linear systems. 

Proof: Let X e Sp{4, R). Then so is X'^ and hence X'^X e Sp{4, R). Let bl(E)l + c(E>i + 
<8) i + e (g) A; be its quaternion representation. Note, as X^X is positive definite, b> 0, while 
c, d, e are pure quaternions. Let i? be a real symmetric, symplectic matrix with non-zero 
trace satisfying H"^ = X^X. Suppose H = al®l+p®i + q^j + r®k. Then equating H"^ 
to X'^X yields the system of equations: 

b = + p.p + q.q + r.r (3.1) 

c = 2ap + 2qxr 
d = 2aq + 2r xp 
e = 2ar + 2p x q 

In addition, as H is symplectic, this system is augmented by the conditions i) a > 0; ii) 
q = iii) 1 = + q.q — p.p — r.r. 

The second equation in the system Equation (3.1), together with 2aq = r x p, yields 
q = j^. So, if we can find a (note, we know that a solution with a > exists, in principle) 
then we obtain q. 

At this point it is convenient to divide the argument into two cases. The case g ^ and 
the case q = 0. Note (as a ^ 0) these are equivalent to the cases d^O and d = respectively. 

The case d ^ 0; To find a, it is first noted that H is symplectic, we have 

1 = a^ + q.q — p.p — r.r 



and hence = + q.q = + This yields a quadratic for a^, viz., 

A 6 + 1 9 d.d 

Note that the discriminant of the above equation is ^ , and this is positive since 

{b+l)"^ — d.d = 6^ + 1 + 26 — d.d can be expressed as 26^ + 26 — c.c— e.e (by using the fact that, 
X'^X is symplectic and hence 6^ + d.d - c.c - e.e = 1). But 26^ + 26 - c.c - e.e = 6^ + d.d + 26 
(since 6^ — c.c — e.e = d.d). But 6^ + d.d + 26 > 0, since 6, being one- fourth the trace of the 
positive definite X^X is positive. Hence both its roots are real. Further, since its coefficients 
change sign (^^ and ^ are positive) both these roots are positive. Pick any one root and 
let a be the positive square root of it. This yields a and hence q. 

To find r and p, one inserts the expression q = into the Equations for c and e and 
uses the vector triple product identities to find: 

2 

c = 2ap + -[{r.r)p - {p.r)r\ (3.3) 
2 

e = 2ar -\ — [{p.p)r — {p.r)p\ 

This would yield a linear system for the unknows p,r in terms of the knowns c and e 
provided we can express p.p.r.r,p.r in terms of e,c. To achieve this, first note that 

c.c = 4a^p.p + A\\ q X r W'^ +8a^p.{q X r) 
e.e = Aa^r.r + AWpx q\f +8a'^r.{px q) 



Hence, 



c.c - e.e = 4a^p.p - Aa^r.r + 4||gxr|p-4||pxg||^ 



Next, using the identity {u x v) .{u x v) = {u.u){v.v) — {u.v){u.v), we find that 

N X r ||2= {q.q){r.r) 



(since q.r = 0) and that 



(since q.p = 0) 



pxq \ f= ip.p){q.q) 



Thus 

c.c — c.e 
p.p - r.r = — — - — 
4a^ — \q.q 

Since, we also know that = o? -\- q.q = \ -\- p.p + r.r, we have that 

6-1 

p.p + r.r = — - — 

So, one gets a Unear system for p.p and r.r, in terms of already determined quantities. To 
find p.r, we compute that 

c.e = {2ap+2q x r).(2ar + 2px q) = 4a^{p.r) + 4(g xr).{px q) 

Hence, 

c.e = 4o^(p.r) + 4[{q.p){r.q) - {q.q){p.r)] = 4(a^ - q.q){p.r) 

Now note that a? — q.q ^ 0. Indeed, if it was zero, then = 3^^. But recall, the quadratic 
for was Wa'^ + d.d = Wa^^. So if = , this would then imply that = ^ , which 
is not true. So, we find p.r = j^^rif^^- This yields p.r. 

Inserting these values for p.p, p.r and r.r into Equations (3.3) yields a linear system for 
the vectors p and r: 



c = ap + (3r (3.4) 
e = pp + 'fr 

with a = 2a + ^=2a+^+ ,ITT ^ > P = , j = 2a+^=2a+^ + 

c.c— e.e 
a(4a^— 4q.q) ' 

The system Equation (3.4) is readily solved. Indeed, the system is invertible as the Cauchy- 
Schwarz inequality reveals that the quantity 017 — is at least 4a^. 

Thus, we have found H. If this H is not positive definite, then the H corresponding to 
the other root of the Equation (3.2) has to be positive definite, since that is the only other 
H which is symplectic, symmetric, and with positive trace satisfying H"^ = X'^X. 

The case d = Q Now Equation (3.2), when d = Q has two roots, viz., (which is stricltly 
positive) and 0. Thus, by picking a = ^^^^ , q = 0, p= ,^ and r = , we find the only 



H which is symplectic, real symmetric, with positive trace and which satisfies H'^ = X^X. 
Thus, by uniqueness, H must be the unique positive definite square root of X'^X. 

In the constructive proof of the above theorem, a key step was an explicit formula for the 
trace of H"^. This, calculation can be used to obtain a complete characterization of positive 
definite symplectic X. Furthermore, it indicates which root of Equation (3.2) to pick to 
compute the desired H. 

Theorem 3.3 Let X = a{l®\)+p^i-\-q®j-\-r^k'be a. real symmetric, symplectic matrix. 
Then it is positive definite iff I) a > 0, II) 2a^ — 2q.q + 1 > 0. In particular, a symmetric and 
symplectic with matrix for which a > 0, g = is positive definite. 

Proof: It is well known that a real symmetric matrix is positive definite iff all its eigenvalues 
are positive. Let, Px{x) = a;" — an-ix^~^ + a„_2a;"~^ + . . . + (— l)"ao be the characteristic 
polynomial of X. Then, by Descartes' rule of signs, for instance, X is positive definite 
iff > for all i = 0,...,n — 1. Since X is symplectic, in addition, its characteristic 
polynomial equals its reverse. Thus, we are guaranteed that Px{x) is of the form Px{x) = 
x'^ — a^x^ + a2X^ — a3X + 1. Now clearly 03 = Tr(X) = 4a. It remains to find 02. Now, 02 = 
i[(Tr (X))2-Tr (X^)]. But Tr {X^) = 4:{a^ + q.q+p.p + r.r). Using, 1 = a? + q.q-p.p-r.r, 
we see that Tr (X^) = Sa^ + S,q.q - 4. This yields, 

Px{x) =x'^ - 4aa;^ + (40^ - 4q.g + 2)x'^ - Aax + 1 

Hence positive definiteness of X is equivalent to a > and 2a^ — 2q.q + 1 > 0. In particular 
if a > and g = 0, X is positive definite. 

Prom this follows a strengthening of Theorem (3.1), which is worth stating separately: 

Theorem 3.4 Let X e Sp{A,R). Then there exist scalars a,wo,t'2, a unit quaternion u and 
pure quaternions p, r satisfying the constraints a > 0, 2a^ — 2q.q-\-l > 0, a? —p.p+q.q — r.r = 
li 9 = and + V2 = 1 such that X admits the following quaternion representation 

X =[u(^{vo + V2j)][a l(g)l+p(g)i + q(^j + r(S)k]. 

Remeirk 3.1 Squeezing Transformations: The authors of [1] argue persuasively that the most 
general squeezing transformation is given by (the metaplectic representation of) a positive- 
definite symplectic matrix. Thus, the above theorem obtains a quaternionic parametrization 
of such transformations for two-mode systems. 



More generally, any application which requires the non-compact part of 5^(4, R) would 
benefit from the parametrization in the above theorem, [1, 7, 8]. Note the matrices which 
are symplectic, but not orthogonal [i.e., the non-compact portion of 5*^(4, i?)] do not form a 
group. Hence, the polar decomposition is arguably the most economical fashion to extract 
the non-compact part of an X G Sp{4:, R). 

Prom Theorem (3.4), it is now clear which root of Equation (3.2) to pick for ensuring a positive 
definite H, when d ^ 0. Clearly, the larger a is, the smaller the q.q is. Thus, the condition 
2o^ — 2q.q + 1 > is more likely to be satisfied by the candidate with the larger a. In fact, 
since there is precisely only one positive definite square root, the H corresponding to the 
larger value of a is positive-definite while the other candidate is not (this can also be directly 
verified). So one should pick the larger value for a, i.e., one should pick a to be the positive 



square root of — — 2* ~ • 

We summarize the above discussion into an algorithm for finding the polar decomposition 
of an X e Sp{4:,R): 

Algorithm 3.1 • 1. Compute directly X^X - this is just the Gram matrix of the 
columns of X. X^X is positive definite (and thus symmetric) and is also symplec- 
tic. 

• 2. Compute the quaternion representation of X^X. This is guaranteed to be of the 
form \ih{l®l) + c®i-\-d®j + e®k (with b € R,c,d,e € P) since X'^X is symmetric. 
Furthermore, these quantities are linear in the entries of X^X. In addition, b, which is 
J the trace of X^X, is positive. 

• 3. Let H he a. matrix which has positive trace, real symmetric and symplectic, and 

which satisfies = X^ X. One such H is the positive definite factor in the polar 
decomposition of X. Let H, which is symmetric, have quaternion representation a(l 
l)+p®i + q®j + r®k (with a G R,p.q.r G P). Steps 4-6 below show how to compute 
a,p, q,r. 

• 4. If = 0, then compute a = ^^^^ , q = 0,p = and r = -y^j- This is guaranteed 



b+l I . d.d 



to be the positive factor in the polar decomposition of X. 



• 5. If 7^ 0, then let a be the positive square root of the larger of the two strictly positive 
roots of the quadratic — + ^ = 0. Define Q = j;^- Find p,r by solving the 
linear system of equations Equation (3.4). 

• 6. This yields H as the unique positive definite square root of X^X and thus the 
symmetric part of the polar decomposition of X. 

• 7. Next compute XH~^. For this no matrix inversion needs to be performed. Instead 
use the fact that = —JiH^J^. This shows that if 



H = 



D 

is the block-form of H, then 

_(JT j^T 

Thus, can be written by inspection. The matrix XH~^ is guaranteed to be sym- 
plectic and orthogonal, and thus will admit the quaternion representation (t;o + i'2j), 
with u a unit quaternion and + ^2 = 1- 



V 



/ 



8. To compute u,vq,V2 use the algorithm described in [2]. 



To give a different perspective on the "paucity" of symplectic, symmetric square roots of 
X^X when 7^ 0, a brief detour will be taken through a calculation which is relevant for 
Theorem (3.5) below as well. 

Thus, suppose that H \s a, real symmetric, symplectic matrix with non-zero trace satisfying 
H'^ = X'^X, and suppose H is another. Then, since H, H are square roots of X'^X [in the 
sense of Definition (2.1)], there exists a real orthogonal matrix U with H = HU . Clearly U 
must be symplectic too. Further the condition, H"^ = reads as HUHU = HH, which, 
since H is invertible, is equivalent to 

UH = HU^ 

Similarly, the condition that H be real symmetric is also equivalent to the same condition 
UH = HU^. Indeed, {HU)^ = HU is equivalent to U^H = HU, which is equivalent to 

H = UHU 



This in turn, is equivalent to UH = HU^ . 

So, considering H as fixed we examine which symplectic orthogonal U lead to UH = HU^ . 
Let U be represented as u®v. Recalling that the imaginary part of v is just V2j, we find 
that UH = a{l (g)l)+/3(8>p + 7(8>i + ^0fc + s(8)l + l(8>t, where the scalar a and the pure 
quaternions /?, 7, (5, s, t are given by the equations 

a = auoVo + (Im u.q)v2 (3.5) 



(3 = uqVqp + wo(Ito u X p) + uoV2r + V2{lm u x r) (3-6) 



7 = uovoq + vo{lm u x q) + av2lm u (3.7) 



6 = —U0V2P — V2{lm u X p) + uovor + vo{lm uxr) (3.- 



s = a(Im u)vq — UQV2q — V2{lm ux q) (3.9) 



t = —[vo{lm u).p + V2{}-'m u).r\i + [auQV2 — VQ^bn u).q\i + [u2(Im u).p — VQ{lm u).r]k (3.10) 

A similar calculation for HU^ reveals that the only non-zero terms in UH — HU^ are the 
following: 

• An 1 i term, viz., —{21m{u).p)vo — {21m{u).r)v2- 

• An 1® j term, viz., 2auoV2 — {2\m{u).q)vQ. 

• An l®k term, viz., 2{lm{u)).pv2 — 2{lm{u).r)vo- 

• An s (8) 1 term (with s G P), viz., 2auolm{u) — 2uQV2q — 2v2\ra{u) x q. 

Equating the first three terms to zero one finds that if vn 7^ 0, then — = = 

Im(tt).r Im(tt).g 

Im(tt).p lao 

This forces Im(zt).r = Im(u).p = ^2 = and \m{u).q = 0. Now as H is symplectic, with 
its q term and its trace is not zero, we have that q is proportional to r x p. So lm{u), being 



perpendicular to p and q, must be parallel to q. But it is also orthogonal to q. So lm{u) = 0. 
Hence u is either +1 or —1. Note that under these conditions the s <S>1 term also vanishes. 
Finally, V2 = forces v to be +1 or —1. So, the only choices for uxv, when vq ^ 0, are 1 (E) 1 
or -1 (g) 1. 

If Vo = 0, then one finds lm{u).p = lm{u).r = 0, again. Thus lm{u) is proportional to 
q. Furthermore, as a ^ and V2 ^ 0, we get Uq = 0. Thus, u is purely imaginary and 
proportional to q. Once again, the s (g) 1 term vanishes uner these conditions. Thus, the only 
choice for u x w, in this case, are ® j or — jj^ J - 

Thus, there are precisely four symplectic, real symmetric square roots of X^X which have 
non-zero trace (two with positive trace, and the remaining two with negative trace). If H 
has positive trace, then it is easy to see (e.g., by inspecting the 1 (8) 1 term in H) that of the 
remaining three candidates only the one corresponding to U = (g) j can also have positive 
trace. 

Calculating the Characteristic Polynomial of a Symplectic X: We now show how to calculate 
the characteristic polynomial of a symplectic matrix, assuming its quaternion representation 
is available, without any need for a determinant calculation. Of course, calculating the quater- 
nion representation is required. But since this provides the polar decomposition also, it seems 
worthwhile to pursue this route for the characteristic polynomial. 

Thus, let X be symplectic, with [u (8) {vq + V2j)][a l^l+p^i + ® j + r ^ k], its 
quaternion representation. This now reads as the expression 

a{liSil)+piSip + 'y<S:j + 5<S:k + siSil + liSit 

The scalar a and the pure quaternions /3, 7, 6, s, t, have exploit expressions given by Equation 

(3.5), . . ., Equation (3.10). 

Since its characteristic polynomial Px{x) equals its reverse, it is of the form Px{x) = 
+ aax^ + a^x^ a^x -h 1. Furthermore, 03 = -Tr(X) and a^. = i((Tr(X))2 - Tr{X'^)). 
Clearly, then 

as = —Aa = —4:{auoVo + (Im u.q)v2) (3-11) 

To find 02 one needs the trace of X"^. A useful observation here is that for this one does 
not need to calculate X"^ . Once X has been found the trace of X"^ is simply 4(a^ + /3./3 + 



7-7 + ^-S — s.s — t.t). 

Now an explicit calculation (which makes repeated use of the fact that u,v are unit 
quaternions) reveals + (3.(3 + 7.7 + 5.5 — s.s — t.t is the expression 

(^o~^2)[9-<i'+(^o~ II ^"tnu \ \^)a^—2{lmu.q)'^]+p.p+r.r—2[{lmu.p)^+{lmu.r)^]+8auoVoV2{lrnu.q) 
Hence the expression for 02 is 



a2 = Sa'^ulvl+Svlilm u.qf+2{v2-vl)[q.q+{ul- \\ Imu | p)a^-2(Im u.g)^]-2(p.p+r.r)+4[(Im ■u.p)^+(Im it.r)^] 

(3.12) 



The above calculations can be summarized in 

Theorem 3.5 Let X e Sp{4:,R) be represented by quaternions as in Theorem (3.4). Then 
its characteristic polynomial is expressible as Px{x) = x'^ + asx^ + a2x'^ + asx + 1, with 03 
and 02 given by Equation (3.11) and Equation (3.12) respectively. 

Remark 3.2 I) In computing the polar decomposition of a symplectic X, one could resort 
to first diagonalizing the symmetric matrix X^X in the form {X^ X)V = D (with V 
orthogonal) and then finding the unique positive definite square root, Y, in the form Y = 
VOiv^. Here is the diagonal matrix whose entries are the positive square roots of 
the entries of the diagonal matrix D. For this, one could use the method of [16], which 
uses quaternion methods for diagonalizing real symmetric 4x4 matrices. Recall that the 
key step (though not the sole step) in this is the following. Given a real symmetric matrix 
a l^l+p'S'i + q'Sij + f^k, one has to compute the singular value decompositon of the 
real 3x3 matrix Z = \p q f], i.e, the matrix whose columns are the vectors in equivalent 
to the pure quaternions p, q, f. In this regard, it is noted that this calculation is somewhat 
simplified for the purpose at hand. Indeed, X G Sp(A,R) implies X^ X € Sp{A,R), i.e., the 
symmetric matrix X^ X is also symplectic, and thus the corresponding terms p, q, r have the 
property that q is orthogonal to p, f. So when computing the SVD of the corresponding Z, 
equivalently when forming the Grammian of the triple of vectors p, q, r, one finds that 



\ 



p.p p.f 



z'^z = 







q.q 







V 



p.r 







r.r 



I 



Hence finding tlie SVD amounts to just diagonalizing tlie 2x2 matrix 

( - -\ 
p.p p.r 

^ p.r r.r ^ 

This can be done by hand. However, partly because, finding the SVD of Z is not the sole 
step, this approach was eschewed in favour of Algorithm (3.1). 

n) Computing the Euler-Cartan Decomposition of 5^(4, R): The Euler-Cartan decomposition 
of Sp{4:,R) asserts that every X e Sp{4:,R) can be factorized as U1DU2, with Ui symplectic 
and orthogonal, for i = 1,2, and D diagonal, positive-definite and symplectic. This factor- 
ization can be explicitly computed starting from the constructive polar decomposition given 
for Sp{4:, R) in this work. Indeed, let X = UP be the polar decomposition. Following the 
technique in part I) of this remark, diagonalize the positive definite symplectic matrix P. It is 
simple and interesting to verify that the diagonalizing matrix V yielded by this procedure is 
actually symplectic and orthogonal. Thus, let V^PV = Dhe the diagonalization of P, with 
V symplectic and orthogonal, and D positive definite symplectic. This yields X = U1DU2 
with Ui — UV, U2 = . Clearly the Ui are symplectic and orthogonal. This then yields (in 
closed form) the Euler-Cartan Decomposition of X. 

4 Conclusions 

This work produced a quatcrnionic representation for the symplectic group 5*^(4, R). Part of 
the utility of this work is that, in attempting to provide explicit formulae for this represen- 
tation, one obtains a very simple and explicit technique for finding the polar decomposition 
of matrices in S'p(4, i?). Since the positive definite factor in this decomposition is one repre- 
sentation of the non-compact part of the symplectic group, this circumstance can be used to 
address applications where this factor is relevant, [1, 7, 8]. It would be an interesting exercise 
to derive other decompositions of S'p(4, R) starting from the representation provided here. 

Finding similar quaternionic representations of matrix groups preserving other bilinear 
forms in dimension four is an interesting question. It remains to see whether these lead 
to as elegant a set of expressions, such as those for the symplectic group. Indeed, it is no 
exaggeration to say that a pivotal role in the results here is played by the fact that the q 
term, in both Theorem (3.1) and Theorem (3.4), is essentially the cross product of the r 



and p terms. Extending such representations to higher dimensions is also open. Quaternionic 
representations are, of course, Hmited to dimension four, just as there is no similar extension of 
the pair of unit quaternions representation to higher dimensional orthogonal groups. However, 
one can hope that in conjunction with either numerical techniques for the symplectic group, 
[4] or methods of Clifford Algebras, [15] these results can be extended to higher dimensions. 
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